Valley Splitting in Si-Inversion Layers at Low Magnetic Fields 



O 

o 

(N 

< 

0\ 



V. M. Pudalov", A. Punnoose^, G. Brunthaler'^, A. Prinz"^, G. Bauer"^ 
" P. N. Lebedev Physics Institute, Moscow, Russia. 
' Weizmann Institute of Sciences, Rehovot, Israel 
Institut fur Halbleiterphysik, Johannes Kepler Universtdt, Linz, Austria 

We report novel manifestation of the valley splitting for the two valley electron system in (100) Si- 
inversion layers at low carrier density. We found that valley splitting causes almost 100% modulation 
of the Shubnikov de Haas oscillations in very low magnetic fields, almost on the bound of the quantum 
interference peak of the negative magnetoresistance. From the interference pattern of oscillations 
we determined the valley splitting in the B = limit which appears to vary within a factor of 1.3 
over the density range (3 — 7) x 10^^ cm~^. Within the same range of densities, level broadenings in 
both electron valleys differ only by < 3%. The latter result shows that the inter- valley scattering is 
not responsible for the strong (six fold) 'metallic-like' changes of the resistivity with temperature. 

PACS: 71.30.-hh, 73.40.Hm, 73.40.Qv 
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The apparent 'metallic-like' temperature dependence 
of the resistivity is in the focus of current research inter- 
est. A number of microscopic mechanisms have been pro- 
posed for explaining the two major features of the resis- 
tivity: (i) strong metallic-like changes with temperature 
|1-S|, and (ii) strong changes with in-plane magnetic field 
1 7 - ni . One of the promising models intensively discussed 
in this connection is based on scattering between different 
sub-bands [^-0. If the carriers mobility strongly de- 
pends on the Fermi energy Ep, the minority subband(s) 
may have a mobility essentially lower than that for the 
majority subband(s). The intersubband scattering may 
thus cause strong changes in the resistivity with temper- 
ature and magnetic field. In addition to the spin-related 
subbands, n— (lOO)-Si system has two minima in the con- 
duction band. The minima originate from six equivalent 
valleys located close to the AT-points in the Brillouin zone. 
Four of them shift up in energy by about (20-40) meV due 
to the confinement potential. At low densities, only the 
two lowest valleys are filled; the sharp Si/Si02 interface 
causes their additional splitting by Ai, > 1 K . There- 
fore, apriory, the remarkable strength of the metallic-like 
conduction in Si-samples might be related to the valley 
multiplicity. 

The 'metallic' temperature changes in the resistivity 
and the 'metal-insulator' transition (MIT) in 2D sys- 
tems 1^,^ are most pronounced in high-mobility samples 
where they take place in the regime of low carrier density 
and strong electron-electron interaction (the ratio of the 
Coulomb interaction energy to Fermi energy, Tg , is of or- 
der of 3-10). The interaction is expected to enhance both, 
spin- and valley-splitting ||l5|; a possibility of a sponta- 
neous polarization, caused by interactions, was discussed 
earlier ||l6|] and was recently recalled again in connection 
with the problem of the MIT in 2D . 

In strong quantizing magnetic fields B± {lUcT ^ 1), 
valley splitting is known to be enhanced by exchange in- 
teraction between energy levels [nsl and was directly mea- 



sured from Shubnikov-de Haas (ShdH), quantum Hall 
(QHE) effect, and chemical potential oscillations. In par- 
ticular, in Ref. ||l^, valley splitting was found to be al- 
most density independent and to increase linearly with 
magnetic field, A„(B) = A° + aB, where a = 0.6K/T 
and the sample dependent A° ?» 2 K. The latter renor- 
malization is intrinsic only to the strong field regime of 
well separated energy levels and is irrelevant to the prob- 
lem of metallic conduction at i? = 0; valley splitting in 
Si-structures at i? ^ remained so far unexplored. 

In the current work we have found that valley splitting 
manifests itself in high mobility Si-samples not only at 
high fields, but also at low magnetic fields, B < 0.8 T, 
giving rise to beatings in the ShdH oscillations. Such 
beatings were discussed theoretically in Ref. [|l5| b but 
have not been observed experimentally. Due to very high 
mobility of the samples we traced the oscillations down 
to the fundamental limit of w 0.2 T, given by the onset 
of the quantum interference (i.e., weak negative magne- 
toresistance) |l^. From this novel beating pattern we 
determined the zero-field valley splitting as a function of 
the carrier density. We found that A„ varies only weakly 
in the range of low densities, n ^ (3 — 7) x 10^^ cm^^. Al- 
most 100% amplitude of modulation in the beating pat- 
tern evidences that the scattering times in the two val- 
leys are about equal. This observation demonstrates that 
the inter-subband scattering mechanism is not responsi- 
ble for the strong 'metallic-like' temperature dependence 
of the conduction in Si-inversion layers. 

We performed measurements on three high mobility 
n— Si-MOS samples selected from three different wafers 
(whose plane coincide with (lOO)-crystal plane to within 
1°): Sill (peak mobility /iP<='^'= = 3.9mVVs at T = 
0.3 K), Sil2 (^P"'*'^ = 3.4) and Sil5 (/xP'^^'^ = 4.0). All 
samples exhibited strong (six-fold) metallic-like fall in the 
resistivity with temperature, and an apparent MIT at a 
'critical density' Hc « 0.9 x lO^^cm"^. 

Figure 1 shows a picture of the oscillations typical for 
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high magnetic fields and for a relatively high carrier den- 
sity, n = 10.24 X 10^^ cm"^. As the magnetic field in- 
creases, ShdH oscillations evolve into the quantum Hall 
effect. In quantizing fields the energy spectrum is 

e^huj,{N+^)±^gyBB±^A4B), (1) 

where g* is the effective g— factor, /i^ is the Bohr mag- 
neton, and > is the Landau level index. 
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FIG. 1. Typical manifestation of valley splitting in di- 
agonal resistivity at high fields and high carrier density. 
Arrows point at the valley gaps u — 5,7, encircled num- 
bers depict filling factors for the cyclotron gaps, and 
V = 6,10 are for the Zeeman gaps. Inset illustrates the 
hierarchy of the energy levels in high fields. 

The inset to Fig. 1 illustrates the energy spectrum 
Eq. (1) in the high field regime. The strongest minima 
in Fig. 1 a (and the largest splittings) at high fields cor- 
respond to the cyclotron gaps Ac = hiVc — g* ^ibB — A^,; 
they occur at integer filling factors = 4, 8, 12..., where 
V = nh/eB±. The series of minima at i/ — 2,6, 10, etc. 
corresponds to Zeeman-gaps, = g* hbB — A^,. Start- 
ing from i? = 6T, valley gaps become resolved and pro- 
duce minima in pxx a^i v — 7, 5, 3, 1. For magnetic fields 
less than 6 Tesla, valley splitting is not seen in Fig. 1 at 
any density, to within 1% resolution. 

Spin splitting (i/ ~ 6, 10, 14, etc) can be seen in Fig. 
1 down to about 1.5 T. For lower fields, the oscillations 
decay exponentially in amplitude, become harmonic and 
obtain the anticipated periodicity ABe/hn corresponding 
to four degenerate levels (2 spins and 2 valleys). 

Unexpectedly, in the low field range, B = 0.2 — 0.3T, 
almost on the border with the quantum interference peak 
in Pxx{B) p9[ , we found beatings in the oscillatory pic- 
ture, demonstrated in Figs. 2. The beatings are ob- 
served in all three samples but only at low densities 



n < 9 X lO^^cm"^. For n < 3 x lO^^cm"^, the de- 
creasing number of oscillations in the interference pattern 
prevented us from making a quantitative analysis. 
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FIG. 2. Shubnikov-de Haas oscillations at low magnetic 
fields and low densities: a)over an extended range of fields, 
for n = 8.39 x lO^^cm"^ ; b) and c) magnified low field 
region of the interference pattern, at two carrier densities, 
n = 8.39 and 7.27 x 10^^ cm~^. Arrows mark the beating 
nodes. 

The interference picture in Figs. 2 corresponds obvi- 
ously to the presence of two oscillatory terms in Pxx{B~^) 
with slightly different periods and almost the same am- 
plitude. We note, firstly, that Zeeman splitting in purely 
perpendicular field in no case may produce such interfer- 
ence. Secondly, since both, cyclotron and Zeeman energy 
vanish as B± decreases, valley splitting should become 
the dominant parameter in low magnetic fields. Indeed, 
based on the measurements of A^, iQ and of 5*-factor 
pO| we estimate 

Ac = hiVc - Az - A^, « 1.8,B±_ - 2, 
Az - g*tJ.BBi_ - A^ « 25^ - 2 

A^«2 + 0.6Si, (2) 

where all energies are in K, B± is in T, Hluc ~ 7 x 
(0.19/m*) « 5K/T [|l5|,|o), and » 2.6K/T |2^] 

for a typical density n = 5 x 10^^ cni^^. 

We presume therefore that the interference is caused 
by valley splitting. With this presumption, we mod- 
eled the oscillatory component of the magnetoresistance 
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SpxxiB±) with Lifshitz-Kosevich (LK) formulae [ piy22| 
for spin-degenerate carriers: 



hn, 



■v-\- 



2eB± 



cos I 2tts——- TTs 

2eB I 



(3) 



Here ny± — (n/2) [1 ± Ay/{2ep)] are the partial popula- 
tions in the ± valleys. The envelope function is 



with 



exp 



2Tr^skT/huj* 



-FAg*m*), (4) 



fi/{2'KTjj) being the quantiun lifc time, 
and T^, the Dingle temperature in each ± val- 
ley. The Zeeman factor, F = cos^nsg* hbB /hujc) — 
cos{nsg*m* /2mQ) ||2l[], and m*{n), tjiq, are the effective 
and free electron masses correspondingly. 



which are at least 3 times lower; however, the oscillations 
vanished earlier. 

We obtained more than a satisfactory fit using Eqs. (3) 
and (4) for modeling over a whole range of densities where 
the beatings are observed; an example is presented in 
Fig. 3. The main result of our fitting, valley splitting 
for two samples is presented in Fig. 4. appears to 
be sample dependent, even for the samples with similar 
peak mobility (same disorder). 

As seen from Fig. 4, over the explored interval of den- 
sities, n = (7.5 — 3) X lO^^cm"^, the changes in A„ are 
about a factor of 1.15 - 1.3, much bigger than the error 
bar, ~ (2 — 5)%. The changes though are not reminis- 
cent at all of a critical behavior, cx; l/(n — no), which 
might be expected for the spontaneous valley polariza- 
tion at a certain density no ric. It is noteworthy that 
over the same interval of densities both, the electron ef- 
fective mass m* and (7*-factor were found ||2^ to change 
monotonically by a factor of about 1.2. 
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FIG. 3. Example of fitting of the measured inter- 
ference pattern with Eq. (2) for the sample Sill, 
n = 6.04 X 10^^ cm"^; T = 0.3 K. 

In modeling, we used two fitting parameters, A„ 
for valley splitting and fq — (t+ -|- t~)/2 for quan- 
tum lifetime. The former parameter determines en- 
tirely the Fourier spectrum of oscillations and nodes 
location, whereas the latter one describes the absolute 
amplitude of oscillations and its monotonic field depen- 
dence. The difference of partial lifetimes for each valley, 
(t^ — T~)/Tq <^ 1 was used for a fine adjustment of the 
modulation depth at the node and is found to be a small 
parameter, < 3%, within the explored range of densities. 
Two other parameters of the energy spectrum, g* (n) and 
m*(n) were determined in Ref. pO| ] vs carrier density. It 
is important that the desired valley splitting is deter- 
mined entirely by the location of the node of oscillations 
on the magnetic field scale (as illustrated in Fig. 3) and 
is almost insensitive to g*, m* and values. 

We emphasize that the node of beating demonstrated 
in Figs. 2 and 3 is the first node; this circumstance en- 
abled us to determine unambiguously the valley splitting. 
The second node might be expected at magnetic fields 
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FIG. 4. Valley splitting versus carrier density for two 
samples. Filled symbols are for Ai,, open symbols are for 
AJJ. The uncertainty is represented either by error bar or 
the symbols size. 

Our measurements were performed in very low though 
finite fields, B « 0.3 T. It is not clear whether the em- 
piric field dependence Eq. (2) remains valid down to 
B = or Ay decreases more rapidly as levels starts 
overlapping in low fields. The accuracy of our data 
was insufficient to determine this experimentally though 
we found the fitting to be better with field dependent 
Ay{B) (as in Eq. (2)). Given Eq. (2) is applicable at 
B — » 0, it contributes about 20% to the obtained A„ 
values. We present therefore in Fig. 4 the results of 
both fitting, with field independent A^ and with field 
dependent A^ — A° -I- 0.6B. The two results have the 
meaning of the upper and lower estimates, correspond- 
ingly, for the true Ay{B — 0) value; their difference is 
about equal to 0.6i?nodc [K/T], where Snodc is the lo- 
cation of the 1st node of oscillations on magnetic field. 
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The sought-for A(_B = 0) value is thus within the inter- 
val from A„ to Aj). Both, A„ and A^ exhibit a weak 
non-monotonic dependence whose origin is unclear. Our 
analysis of the ShdH oscillations in the two- valley system 
performed up to = 5 didn't reveal any deviation from 
the LK-formulae. This justifies the assumption of the 
conventional Fcrmi-liquid behavior which we used in the 
analysis above. 

An additional important result follows from our fit- 
ting. We found the quantum lifetime Tq to be almost 
the same in the two electron valleys, the maximal dif- 
ference being only (t+ — t~)/t^ « 3% for the density 
of 3 X lO-'^'^cm"^. The proximity of values mani- 
fests in Figs. 2 as deep falls in the oscillations ampli- 
tude at the beating nodes. The transport time (momen- 
tum relaxation) for Si-inversion layers is very close to 
Tq H^]; it is therefore very likely that the partial trans- 
port times and partial mobilities in both valleys are al- 
most equivalent. This is in contrast to the big, six-fold, 
changes of the resistivity with temperature at zero field, 
(p(Thigh) — p(7iow)) / p{T) which the sample exhibit at 
the same density p^ ]. From the proximity of the quan- 
tum life times in both valleys, we conclude therefore that 
the inter-valley scattering mechanism | [T^ is not respon- 
sible for the strong metallic-like temperature variation of 
the resistivity in Si-inversion layers. 

To summarize, in high mobility (lOO)Si-inversion lay- 
ers, the system which exhibits very strong 'metallic-like' 
features in conduction, we observed a novel manifesta- 
tion of valley splitting: it causes unexpected beatings 
in Shubnikov-de Haas oscillations in low magnetic fields, 
(0.15 - 0.4) T, right on the bound of the negative mag- 
netoresistance peak. From the beatings pattern of os- 
cillations we determined the valley splitting A„ in the 
_B = limit. We found that A„ varies with density 
rather weakly and doesn't display a critical behaviour 
in the range of densities (3 — 7.5) x lO^^cm"^ or = 
3-5. We determined also the individual quantum life- 
times, in both valleys, which appear to differ by less 
than 3%. This insignificant difference in Tq demonstrates 
that the semiclassical mechanism of mobility changes re- 
lated to the inter-valley scattering is not the origin of 
the strong metallic-like temperature dependence of the 
resistivity in Si-inversion layers. The findings set novel 
constraints on the microscopic models developed to ex- 
plain the 'metallic'-conduction. 
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I. APPENDIX: ENERGY SPECTRUM IN 
SI-INVERSION LAYER AT LOW MAGNETIC 
FIELDS 

The results on the valley splitting presented above shed 
a light on the energy spectrum of the two-valley elec- 
tron system in (100) Si- inversion layers in low fields. We 
plotted in Fig. 5 the energy for three lowest Landau lev- 
els, = 0, 1, 2, according to Eqs. (1) and (2) (from the 
main section). In high fields, the cyclotron splitting is 
the largest and the sequence of energy levels corresponds 
to the inset to Fig. 1. Counting from the lowest Fermi 
energy, the 1st gap in the energy spectrum (filling factor 
v = 1) corresponds to valley splitting; v = 2 is for the 
reduced Zeeman splitting, = 3 is again for the valley 
splitting and z/ = 4 is the reduced cyclotron gap. This 
picture was verified in numerous experiments ||l|,^. 

Due to the nonzero valley splitting at i3 = 0, as mag- 
netic fields decreases, the energy levels start crossing each 
other at B ^ 1 T . In the region of crossing, this single- 
electron picture fails and the repulsion between levels 
should be taken into account. There is almost no ex- 
perimental data available on the details of the energy 
spectrum in Si in low fields. It is well known only |^|^] 
that in low fields/low density regime quantum oscilla- 
tions are missing for v = 3, 4, 5 and remain pronounced 
for = 1,2 and ior i' — 6 and 10; this result was recently 
confirmed in ref. [f7|. 
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FIG. 5. Schematic magnetic field dependence of the en- 
ergy spectrum (in the single-electron approximation) for 
two valley electron system in Si and for three lowest Lan- 
dau levels. The parameters of the spectrum correspond 
to Eqs. (1) and (2); = 2K, as for the sample Sil2. 



to Fig. 5) a number of energy levels coincide. We don't 
think though that the simplified single-electron picture 
of Fig. 5 may predict exact location of the level crossing 
regions. 

The identification of each energy splitting in low 
field/low density regime is not transparent and requires 
more thorough theoretical calculations of the spectrum 
and much more detailed experiments. As an example, 
on the n — B plane in Fig. 6 we reproduce from Ref. |^ 
a Landau fan diagram of oscillations in the low density 
limit; similar data were reported recently in Ref. 0|. In 
Fig. 6, the y02;2;-minima a,i v = Q can be traced with- 
out interruption down to the field B = 0.55r, through 
the 'metal-insulator' boundary, however the origin of this 
splitting may change from 'Zeeman gap' at high fields to 
any other one, e.g. 'valley gap' or 'cyclotron gap' at low 
fields. By now, there is no firm physical background to 
identify v — Q oscillation in low fields with spin split- 
ting. For example, one can see from the inset to Fig. 5, 
that in the single-electron approximation the v — Q gap 
&i B = 0.6 T corresponds to the transition with chang- 
ing valley (but is reduced by the spin flip energy) and 
V = A gap is the combined gap with changing both, the 
cyclotron number and the spin projection. 




FIG. 6. Landau fan diagram for the four different 
Piijj-minima, reproduced from Fig. 2 of Ref. |^ . Full dots 
show MIT-boundary, defined according to vanishing ac- 
tivation energy. Other symbols designate location of the 
Piijj-minima on the n — B-plane. 



In fields B < 0.7 T, the sequence of energy levels 
changes as shown in the inset to Fig. 5. For example, 
&t B = 0.6 T, the 1st splitting, i/ = 1, corresponds to 
the Zeeman gap, v — 2 la the cyclotron gap (reduced by 
the spin flip energy), z^ = 3 is a combination of the val- 
ley and cyclotron splitting etc. In this simplified picture, 
at some field values [B = 0.45 T and 0.25 T in the inset 



In Ref. 0, from the estimated ratio of the two gaps, at 
V — & (which authors of Ref. |^] assumed to be the spin- 
gap) and at v = A (was assumed to be the cyclotron gap), 
a conclusion was drawn on the unexpectedly strong en- 
hancement of the (7*-factor at low electron densities. Tak- 
ing into account vanishing Ac in the crossing region and 
ill-defined splittings at low field, we find this argument 
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unjustified. Recent direct measurements of the Zeeman 
energy in low fields vs carrier density |^ show that the 
Zeeman energy increases smoothly as density decreases 
down to Tie, in agreement with the anticipated Fermi liq- 
uid renormalization and with recent results by Okamoto 
et al. without any unexpected deviations or diver- 
gence in the range of densities n = (1 — 90) x lO^^cm^^. 

It is noteworthy that for the weak field interference 
pattern which we analyzed above, the interpretation of 
each oscillation is not important and docs not influence 
on the conclusion about the valley origin of the interfer- 
ence and on the results of our analysis. 
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